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ABSTRACT 

We investigate the [On] emission line properties of 18,508 quasars at z < 1.6 
drawn from the Sloan Digital Sky Survey (SDSS) quasar sample. The quasar sample 
has been separated into 1,692 radio-loud and 16,816 radio-quiet quasars (RLQs and 
RQQs hereafter) matched in both redshift and i'-band absolute magnitude. 
*^ We use the [On]A3726 + 3729 line as an indicator of star formation. Based on 

these measurements we find evidence that star-formation activity is higher in the 
RLQ population. The mean equivalent widths (EW) for [Oil] are EW([On]) RL = 
7.80 ± 0.30 A and EW([On]) RQ = 4.77 ± 0.06 A for the RLQ and RQQ samples 
respectively. The mean [On] luminosities are log 10 [L([Oii]) R L/W] = 34.31 ± 0.01 and 
log 10 [L([Oli]) RQ /W] = 34. 192 ±0.004 for the samples of RLQs and RQQs respectively. 
Finally, to overcome possible biases in the EW measurements due to the continuum 
emission below the [On] line being contaminated by young stars in the host galaxy, we 
use the ratio of the [On] luminosity to rest-frame i'-band luminosity, in this case, we 
find for the RLQs log 10 [L([On]) RL /L opt ] = -3.89±0.01 and log 10 [Z([Oll]) R Q/Z opt ] = 
—4.011 ± 0.004 for RQQs. However the results depend upon the optical luminosity of 
the quasar. RLQs and RQQs with the same high optical luminosity log 10 (L op t/W) > 
38.6, tend to have the same level of [On] emission. On the other hand, at lower optical 
luminosities rog 10 (L opt /W) < 38.6, there is a clear [On] emission excess for the RLQs. 
As an additional check of our results we use the [Oin] emission line as a tracer of the 
bolometric accretion luminosity, instead of the z'-band absolute magnitude, and we 
obtain similar results. 

Radio jets appear to be the main reason for the [On] emission excess in the case of 
RLQs. In contrast, we suggest AGN feedback ensures that the two populations acquire 
the same [On] emission at higher optical luminosities. 
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1 INTRODUCTION 

1.1 The link between accretion activity and 
star-formation activity 

It is now clear that black-hole growth, and associated active 
galactic nuclei ( AGN), and the growth of the host galaxy are 
conn e cted (e.g. Magorrian et al.| 19981: Ferrarese fe MerrittJ 



l200d : iGebhardt et al.ll200d ; lHaring fc Rbdl2004l h Further 
more, simulations assume that the triggering of the 
main phase of AGN activity in gas-rich mergers will 



be accompanied by a major galaxy - wide starburst (e.g 



Pi Matteo et al.l 120051; ISpringel et al.1 120051 : iHopkins etall 
ajj 



20081 : ISomerville et al.ll20oi ) 



Email: c.kalfountzou@hcrts.ac.uk 



In order to understand the link between the two 
processes and assess the possibility of the two occur- 
ring concomitantly, it is important to quantify and con- 
strain the star-formation activity in quasar host galaxies. 
From early studies on the star formation - quasar con- 
nection, it was rea lized that the evolving l uminosity den- 
sity o f quasars (e.g iBovle fc Terlevicbi ri998: Richards et al.l 
120051 : [C room et al.1 1200911 and the cosmic star-formation 
history (e.g. iMadau et al.lll99rj; IHopkins fe Beacoml l2006t) 
are similar (e.g. Franceschini et al. 19991 : Archibald et al.l 
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l2002h . Supporti ng evidence has also r eported the pres- 



ence of dust dArchibald et al 
Reuland et al.l 12004; Stevens et al 



200ll; iPage et all l200ll ; 
' 20051). cold gas (e.| 



Wilson fc Colbertl Il995l ; ISikora et aU 120071 ; but see also 
Garofalo et alj 120100 ■ Apart from the central black hole, 



Evans et all l200ll; IScoville et al.l 120031 ; IWalter et all 12004 



Emonts et alj 201l|) and young stars (e.g. Tadhunter et alj 



2005l ; lBaldi fc Capettill2008l ; lHerbert et alj|2010h in powerful 
radio-loud AGN, which does not reconcile with the picture 
of low star-forming ac tivity in the ellipti cals that usually 
host these objects (e.g. iDunlop et alj|2003h . 

Many studies have also attempted to determine the 
star-formation a ctivity in quasar ho st galaxies using opti- 
cal colours (e.g. Sanchez et al. 20041') or spe c troscopy (e.g. 
Trichas et ail 120101 : iKalfountzou et alj 12013 ; iTrichas et all 



20121 ). The star formation and quasar connection has also 
been investigated in recent work (e.g. Trichas et al.l 120091 ; 



ISerieant et aTll2010l ; iBonfield et aHl201ll ) using far-infrared 



and submillimiter data from the Hers chel Astrophysical Ter- 
ah ertz Large Area Survey (H-AT LAS; lEales et afll20ld ) and 



by iHatziminaoglou et all (|2010h using the Herschel Multi- 
Tiered Extragalactic Survey CHerMES; [Oliver et al1l20ld ). 
showing that the far-infrared emission is weakly correlated 
with the accretion rate for optically s electe d quasars. In com- 
plementary work, lHardcastle et all (|2010l ) find that the far- 
infrared emission in radio galaxies appears consistent with 
the population of non-active galax i es, at least at low radio 
luminosities, while ISevmour et all (|201ll ) suggest that the 
star-formation activity in radio galaxy hosts increases with 
redshift. 



1.2 The radio- loud — radio-quiet dichotomy 

The radio-loudness dichotomy constitutes a subject that has 
attracted considerable attention over the years. RLQs are of- 
ten defined to be the subset of quasars wit h a radio-loudness 
satisfying Ri > 10, where Ri = LI - 5GHz ) ||Kellermann et ail 

i(4000A) 

1989) is the ratio of monochromatic luminosities measured 
at (rest frame) 5 GHz and 4000 A. The RQQs must mini- 
mally satisfy Ri ^ 10. It has long been claimed that RLQs 
are only a small fraction ~10 per cent of all quasars (e.g. 
Ilvezic et alll2002l ). with this fraction possibly v arying with 
both luminosity and redshift (I Jiang et alj|2007l ). 

One of the main difficulties in studies aimed at com- 
paring the properties of the two populations is the low 
fraction of radio- loud sources, which leads to difficulty in 
defining statistically meaningful samples of RLQs. With the 
advent of the SDSS database, coupled with the Faint Im- 
ages of the Radio Sky at Twenty cm survey (FIRST survey; 
iBecker et al1ll995l ). it became possible to select a large and 
more complete sam ple of radio- l oud AGN. Recent SPSS- 
based studies (e. g. Ilvezic et all |2002| ; IWhite et all 120071 ; 
lLabita et alj|2008l ) confirm the presence of a dichotomy in 
the radio properties of quasar populations. 

While the definitive physical explanation of this di- 
chotomy remains confused, a large number of models 
have been put forward to explain it. Some of them 
are based on the idea that RLQs re quire more massive 



central black holes th an R QQs (e.g . Dunlop et all 120031 ; 
iJarvis fc McLurel 2002; also IShankarll2010l find this to be 
redshift dependent) and it has also been suggested that 
RLQ s host more rapidly s pinni n g black holes than RQQs 
(e.g. iBlandford fc Znaiekl Il977l ; iPunslv fc Coronitil Il990l : 



differences between the two populations have also been 
found in th e ir en vironment and host galaxies. Recently, 
iFalder et all (|201Gh showed that radio-loud AGN appear to 
be found in denser environments than their radio-quiet coun- 
terparts at z ~ 1, in contrast with prev ious studies at lower 
redshifts fe.g. lMcLure fc Dunlod feoOl). However the differ- 
ences are not large and may be partly explained by an en- 
hancement in the radio emission due to the confinement of 
the r adio jet in a dense environment (e.g. lBarthel fc Arnaudl 
Il996t ). 

At lower redshift, the SDSS has also been used to in- 
vestigate differences in the host galaxies, star formation ac- 
tivity and e nvironments for well-defined samples of radio 
galaxies (e.g.lKauffmann et alll20oA 12004 iBest et alj|2005l ; 
iBest fc Heckmanll2012T)7 Of particular r elevance to the work 
presented here. iKauffmann et all l|2008h suggest causal con- 
nections between the ability of an AGN to produce radio 
emission and the environment in which is resides, with ra- 
dio AGN being found in denser regions. This is in qualitative 
agreement with work on higher redshift quasars discussed 
above. However, we note that th e sample radio galaxies anal- 
ysed bv lKauffmann et all (|2008l ) generally lie at much lower 
radio luminosities than the q uasars under inve s tigatio n in 
this paper, and as noted by K auffmann et all (|2008h the 
radio luminosity appears to be loosely correlated with the 
accretion rate and as such it is difficult to quantitatively 
compare this work with the studies of higher redshift, higher 
luminosity, quasars. 

If the radio-loudness is due to the physics of the cen- 
tral engine and how it is fueled, and the environment plays 
a relatively minor role, the quasar properties may be con- 
nected with the st ar formation in their host galaxies (e.g. 
iHerbert et"aHl201Ch . On the one hand, AGN feedback could 
be stronger in the case of the RLQs due to their higher 
black hole masses and therefore potentially higher radiation 
field, reducing the star-formation rate compared to RQQs 
while on the other hand radio jets could increase the star- 
form ation activity compressing the inter galactic medium 
(e.g. ICroft et al.|[2006l ; ISilk fc Nusserll2010l ). 

In this paper, we present a sample of 18,508 quasars, 
spectroscopically confirmed by the SDSS and with radio 
data drawn from the FIRST survey. We use the quasar 
sample to investigate any possible differences in the star- 
formation activity of RLQs and RQQs using [On] emis- 
sion as a star formation tracer. In Section 2 we describe 
our sample, our selection criteria and our measurements. 
Section 3 presents our results and in Section 4 we discuss 
plausible explanations. Finally, in Section 5 we provide a 
summary of our main results. Throughout this paper, we 
adopt the following values for the cosmological parameters: 
H = 70 km s _1 Mpc -1 , fl M = 0.3 and D. x = 0.7. 



2 THE SAMPLE 

Our quasar sample is drawn from the SDSS Quasar Cat - 
alogue seventh data release (DR7) (jSchneider et al.ll2010l ). 
It contains a total of 105,783 spectroscopically confirmed 
quasars brighter than Mi — —22.0 with at least one emission 
line with full-width half maximum (FWHM) greater than 
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1000 km s _1 or with complex absorption leatures. The full 
range of redshifts is z — 0.065 to z = 5.46, with a mean value 
of z ~ 1.49. For the aims of this paper, where we investigate 
whether there are any differences in the [Oll]A3726, 3729A 
emission line doublet between radio-loud and radio-quiet 
quasars, we restrict the redshift range to z < 1.6 and im- 
pose a signal-to-noise ratio in the equivalent width of the 
[Oil] emission line which must be greater than 3. 

In order to measure the radio emission we cross match 
the SDSS quasars with the Fa int Images of the R adio Sky 
at Twenty cm survey (FIRST; iBecker et al.lll995T ). Sources 
are extracted usin g an elliptical Gaussian fitting procedure 
l|White et al.lll997r i with a 5cr detection limit of ~ 1 mjy. To 
match the quasars to the FIRST radio catalogue, we used 
a 5" arcsec radius from the optical counterpart. However, 
as quasars often show multiple radio components we cre- 
ated a second subset found using a 30" match radius where 
the sub-structured sources were collapsed into single objects 
having radio fluxes equal to the sum of the flux of the var- 
ious components and the position of the final single-source 
being specified by SDSS. With this method we found 9,133 
quasars with matching radius less than 5" and 321 extended 
quasars. However, limiting the sample to only the radio de- 
tected sources, we ignore the main bulk of the population. 
Therefore, to study the entire quasar population in the radio, 
we also stacked the radio images and detected the median 
flux densities of un detected sources, following the method of 
IWhite et al.l (|2007l ). 

We separated the quasar sample into RLQs and RQQs 
based on the radio-loudness, Ri, parameter. The traditional 
radio-loud/quiet division is taken to be the radio-to - optica l 
flux ratio of 10. In this work, following llvezic et al.l l|2002h . 
we compute the radio-to-optical flux ratio, Ri, using the 
radio flux density at 1.4 GHz and the i'-band (A; = 7480 A) 
flux density from the SDSS. 

In order to check that the radio maps from the FIRST 
survey do not miss a significant fraction of extended emission 
around the quasars in our sample due to the relative lack 
of short baselines of the VLA in B-array, we also estimate 
the Rj parameter using NRAO VLA Sky Survey (NVSS; 
ICondon et ailll998h . We find no significant difference in the 
-Ri.4GHz between FIRST and NVSS, with the distribution 
remaining statistically unchanged. This suggests that the 
radio emission from the quasars is generally fairly compact. 

We must also ensure that the two subsets of RLQs and 
RQQs have the same optical luminosity and redshift distri- 
bution, so that any differences found in the [On] emission- 
line properties for the RLQ and RQQ samples are not due to 
differences in their optical properties. As redder passbands 
measure flux from the part of the spectrum relatively in- 
sensitive to recent star formation and also suffer less dust 
extinction, we use the i'-band magnitude to estimate the op- 
tical luminosity. We therefore randomly remove RQQs from 
our parent sample in order to force the same distribution in 
redshift and i'-band luminosity for both the RQQ and RLQ 
subsamples. It is well known that the continuum luminos- 
ity in quasars is a good tracer of the accretion luminosity, 
as t he emission is coming directly from the accretion disk 
fe.g. lReeJT984h . As a result, our matched samples of RLQs 
and RQQs in optical luminosity and redshift should be also 
matched in accretion luminosity. 

Figure [T] shows the final distribution of the 1,692 RLQs 
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Figure 1. Distribution of RLQs (red) and RQQs (grey) in lumi- 
nosity and redshift space. The left and bottom panels show the 
normalized optical luminosity and redshift histograms. 



and 16,816 RQQs in the redshift - optical luminosity (L op t — 
z) plane and their respective histograms. A two-dimensional 
Kolmogorov-Smirnov (K-S) test applied to the L opt — z dis- 
tribution shown in Figure [1] shows that the distribution of 
the RLQs and RQQs samples are indistinguishable with 
a probability that they are drawn from the same under- 
lying distribution p = 0.99. An one-dimensional K-S test 
applied to the log 10 L op t distribution gives: D = 0.01 and 
p — 0.99 while the one-dimensional K-S test to the red- 
shift distribution gives: D — 0.29 and p — 0.97. The mean 
redshift for the two populations is z = 0.83 for the RLQs 
and z = 0.84 for the RQQs, while the mean luminosity 



log 10 (L op t/W) = 38.30 for both. 



2.1 Spectral measurements 

In addition to the [On] measurements, which are used to 
investigate any possible difference in the star formation in 
RLQs and RQQs, we are also interested in the H/3 and Mgn 
emission lines and the continuum luminosity at 3000 and 
5100 A because they hav e been calibrated as vi rial black 
hole mass estimators (e.g. iMcLure fc Jarvisll2002l ) and the 
[Oiii] emission line which has been use d to trace the a ccre- 
tion luminosity in previous studies (e.g. ISimpson|[2005T ). 

Based on the FWHM measurements we separate the 
emission lines into narrow (FWHM < 1000 kms" 1 ) and 
broad (FWHM > 1000 km s -1 ) lines. We fit the narrow 
emission lines with a single Gaussian while broad lines are 
fitted with two/three Gaussians. The noise in the spectra 
is measured using emission line free regions, which cover up 
to 20 A, around each of the emission lines we are interested 
in. The continuum level was removed by fitting power-law 
continua /„ oc v a , where a = const. However, even in the 
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RA=228.Q74, DEC= 51 .783, log 10 (L op /W)=37.33 
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Figure 2. Example spectra (insets) drawn from our samples of RLQs (upper panels) and RQQs (lower panels) with their [On] emission 
line (main panels). The optical luminosity of the objects is shown above each plot with the luminosity increasing from left to right. The 
vertical line in the inset identifies the [On] emission line. The dashed black line in the main panels shows the fitted continuum level and 
Gaussian fit to the emission line. The noise characteristics are also shown for completeness. 



optical-UV range (1200 - 10000 A), the continuum spectrum 
can not always be represented by a single power law. For this 
reason we estimate the local power-law shape of each quasar 
in four ranges, namely, the rest-frame intervals, 1400 - 2200, 
2150 - 3200, 2950 - 4300 and 3900 - 5500 A dNatali et at J 
1998). Figure [2] shows some examples of the spectra and the 
[On] emission line fits for RLQs and RQQs with different 
optical luminosities. 

As a check of our automated method, we compare the 
results of our automatic fitting code with manual measure- 
ments. To perform the test we chose a subsample of 100 
sources with available Mgn, H/3 and [On] emission line and 
find that the equivalent widths are indistinguishable with a 
root-mean square scatter of 0.18, 0.09 and 0.18 around the 
linear 1-1 relation for the Mgn, H/3 and [On] emission lines 
respectively. 



3 RESULTS 

In this section we investigate whether there are any de- 
tectable differences in the [On] emission line properties of 
the RLQs and the matched sample of RQQs, using the dis- 
tribution of the [On] equivalent width, the [On] luminosity 
and the ratio of the [On] luminosity to the accretion lumi- 
nosity (L([Oii])/L op t). 



3.1 The [Oil] emission-line distribution 

Figure [3] presents a comparison of the [On] emission line 
properties of the RLQs sample (red histograms) and the 
matched sample of RQQs (black histograms). The [On] 
emission line properties that are compared include the distri- 
bution of [On] equivalent width, [On] luminosity and [On]- 
optical luminosity ratio (L([Oii])/L op t), where L opt is deter- 
mined using the i'-band. 

The mean equivalent widths are EW([Oii])r,l = 



7.80 ± 0.30 A and EW([Oii])rq = 4.77 ± 0.06 A for 
the RLQ and RQQ samples respectively. The mean 
[On] luminosities are log 10 [L([Oii])r L /W] = 34.31 ± 
0.01 and log 10 [L([On]) RQ /W] = 34.192 ± 0.004 for 
the samples of RLQs and RQQs respectively. Finally, 
in the case of the L([Oii])/L pt ratio, we find for 
the RLQs log 10 [L([Oll])/L op t]RL = -3.89 ± 0.01 and 



l°gio[£([Oii])/£o P t]itQ 



-4.011 ± 0.004 for RQQs. We 



therefore find that the RLQs have different mean values 
for RQQs in all of the parameters. Recall that the RLQ 
and RQQ samples are matched in L op t, so these results im- 
ply that RLQs have higher [On] emission for a fixed value 
of I/opt- We have applied an one-dimensional K-S test for 
each measurement related to [On] for the RLQs and RQQs. 
The K-S test returns a probability of p = 1.12 x 10" 27 , 
p = 3.24 x 10~ 22 and p = 8.77 x 10" 23 and a test statis- 
tic of D = 0.15, D = 0.13 and D = 0.14 for EW([Oii]), 
log 10 [L([On])] and log 10 [Z/([Oll])/Z/ O pt] respectively; thus we 
can reject the null hypothesis that the various measurements 
related to the [On] emission are drawn from the same un- 
derlying distributions. Based on the distribution histograms, 
the RLQs have higher [On] emission than that of their RQQ 
counterparts. 



3.2 The relation between the [On] emission line 

and Li.4ghz 

Previous studies have found a connection between the nar- 
row emission lin e luminosity and the radio luminosity in ra- 
dio g al axies (e.g.lRawlings et a l. 1989; Rawlings fc Saunders 



1991 ; iMcCarthvl 1 19931 ; IWillott et all 1 19991 ; Ijarvis et al 



20011 ). The general assumption here is that the narrow line 
region (NLR) is predominantly photoionised by the quasar 
and the narrow line luminosity is telling us about the intrin- 
sic quasar luminosity. 

In Figure [4] we show the [On] emission line luminosity 
against 1.4 GHz radio luminosity. The strong correla- 
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Figure 3. Comparison of the log 10 [EW([On])], log 10 [L([On])] 
and log 1 g[L([On])/L pt] distribution for the sample of RLQs 
(red) and the matched sample of RQQs (black). The dashed lines 
show the mean values for each population. 



tion between narrow emission line luminosity and radio 
luminosity is clearly seen. The trend fit to the whole 
quasar sample us ing ordinary least squares (OLS) bisector 
l|lsobe et al.lll990l ') gives a correlation of log 10 [L([Oll])/W] = 
(18.94 ± 0.15) + (0.67 ± 0.01) Iog 10 (Li.4GH*/W Hz" 1 sr" 1 ). 
Splitting the sample into the RLQ and RQQ sub-samples 
the equations for the two populations are: 

log 10 [L([On])] = (19.20±0.47) + (0.61±0.02) log 10 (Li. 4 GHz) 

and 

log 10 [L([On])] = (14.88±0.14) + (0.85±0.01) log 10 (Ia.4GHs) 



for RLQs and RQQs respectively. 
The positive correlation of 



[On] emission line and 



radi o luminosity is in agreement with previous studies 
(e.g.lWillott et al.lll999l ; Ijarvis et al.ll200ll ) . IFernandes et all 
performed an analogous study for powerful radio 
galaxies at 151 MHz and they found a strong positive re- 



34 35 
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Figure 4. Luminosity at 1.4 GHz (Li.4GHz) versus luminosity 
of [On] [L([On])] for the RLQ and RQQ subsamples. Red dots 
represent the RLQs and grey the RQQs. Black stars show the 
mean values for each optical bin for RLQs while black squares do 
the same for RQQs. 



lation L([0n]) oc L^*™ 1 . The slightly different slope 
may be explained by the di fferent radio fre q uenci es that 
are used - 151 MHz for the IFernandes et al.1 (120111 ) study 
and 1.4 GHz in our case - or due to the different selected 
source s - powerful radio galaxies for the IFernandes et al.l 
l|201ll ) study and quasars in our study. However, it is clear 
that the relation between radio luminosity and [On] lumi- 
nosity is different for the RQQs compared to the RLQs; this 
again leads to the suggestion that an additional effect is at 
work for ionizing [On] at least. 



3.3 The relation between the i([On]) emission 
line and L opt 

In Figure [S] we show the [On] emission line luminosity 
versus the optical luminosity. A strong correlation between 
the narrow emission line luminosity and the optical lumi- 
nosity is clearly seen with the OLS bisector fit for the two 
populations being: 

log 10 (L opt ) = (3.61 ± 0.64) + (1.01 ± 0.02) log 10 [L([On]) RL ] 
and 

log 10 (L O pt) = (5.76 ± 0.23) + (0.95 ± 0.01) log 10 [L([On]) RQ ] 

for RLQs and RQQs respectively. It is worth noting 
that at high optical luminosities [log 10 (L O pt/W) > 38.6] 
the RLQs and RQQs show a similar distribution in 
log 10 [L([Oll])]. However, at lower optical luminosities, the 
RLQs typically have greater [On] luminosity than the 
RQQs. 

For the quasars in our sample, the black-hole mass 
is strongly correlated with the optical luminosity. Us- 

1 For quasars up to redshift of 1.6 we have used the 
iMcLure fc Jarvisl 120021 1 relations based on the FWHM of H/3 
and Mgll emission lines to estimate the black hole masses. 
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Figure 5. Optical luminosity, L op t, versus [On] luminosity, 
L([On]). Symbols are as in Figure[4] 

ing the L opt — Mbh correlation for the entire sample (Mbh oc 
L ij.8±o.oa^ Lopt = 10 38.6 w corresponds to Mbh = 1Q 8 Mq . 

As suggested by previous studies RLQs tend to have black 
holes more massive than 10 8 Mq. Above the 10 38 6 W opti- 
cal luminosity, all the RQQs in our sample have black holes 
more massive than 10 s Mq , although RLQs continue to have 
a higher mean black-hole mass. 

Figure [6] provides a clearer view of all these parameters 
where we show the optical luminosity versus [On] equiva- 
lent width for the RLQs and RQQs, with black-hole mass 
indicated using a colour scale. It is obvious that the two sam- 
ples tend to have the same equivalent width distribution for 
log 10 (L op t/ W) > 38.6, whereas at log 10 (L opt / W) < 38.6 
there appears to be a pronounced difference. A possible ex- 
planation for these differences between the two populations, 
and for different optical luminosities (or black hole masses), 
could be two different mechanisms that are correlated with 
star formation. At L op t = 10 38 W, the [On] emission in the 
RLQs is roughly 0.4 dex higher than that in the RQQs, 
corresponding to a star-formation rate of ~ 4 Mq yr _1 
(|Gilbank et al.ll2010T ). To explain the [On] emission for op- 
tical luminosities lower than 10 38,6 W we need a positive 
mechanism that increases the star formation in RLQs but 
not in RQQs. The obvious candidate mechanism is that the 
jets observed in RLQs but not in RQQs play a critical role. 
On the other hand, for the higher optical luminosities, we 
need a mechanism that decreases the [On] emission in RLQs 
so that the excess with respect to RQQs is destroyed. As 
the optical limit is correlated with the black hole mass, the 
explanation we need could be provided by negative AGN 
feedback. Both of these mechanisms and their connection 
with [On] emission are discussed in Section 5. 

3.4 [Oni] as a tracer of accretion luminosity 

In this work we use the i'-band magnitude in order to es- 
timate the optical luminosity. We have selected the i'-band 
as an AGN tracer due to the fact that it lies at the red- 
der part of the optical spectrum which is less sensitive to 
the recent star formation. However, any enhancement in the 
[On] luminosity needs to be separated from any possible gen- 
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Figure 6. Optical luminosity, logio(L op t), versus the logarithm 
of the [On] equivalent width, log 10 [EW([Oll])]. The colour scale 
corresponds to the black hole mass estimate. The stars represent 
the RLQs and the squares the RQQs. 

eral enhancement in the emission-line lum inosity in radio- 
loud quasars (e.g. iKauffmann et al1l2008h . In this section, 
we therefore test whether our results are robust to another 
tracer of the accretion luminosity, namely the [Om]A5007 
emission line. The [Om]A5007 line is usually the strongest 
emission-line in the optical spectra of type-2 AGN, and it is 
believed to be ionized predominantly by the photons ema- 
nating from the accretion disk. The [Oni] emission line was 
also used bv lSilverman et all i|2009l) to statistically remove 
the component of the detected [On] emission line that can 
be attributed to the AGN. In order to test the robustness 
of our results, we use the [Oni] luminosity instead of the 
optical i'-band lumin osity and conduct a similar analysis. 
However, we note that iFine et al l l|201lf ) demonstrated that 
[Oni] may not be a robust tracer of the AGN bolometric 
luminosity and may be affected by orientation. 

In our sample of 18,508 QSOs, 1,629 (1,364 RQQs and 
265 RLQs) have the [Om] emission line accessible in the 
optical spectra. In order to ensure that the [On] emission 
excess is not biased due to the selection of optical luminos- 
ity as an AGN tracer, we follow the same method of the 
matched subsamples using the [Om] luminosity as a tracer. 
We separate the sample with available [Om] emission line 
measurements into RLQs and RQQs with indistinguishable 
[Om] and optical luminosity distributions. The subsamples 
contain 237 RLQs and the same number of RQQs. The 
mean values of the [Om] luminosity for the two popula- 
tions are < log 10 [L([Olll])R L /W] >= 34.957 ± 0.011 and 
< log 10 [L([Olll]) RQ /W] >= 34.958 ± 003, while the one- 
dimensional K-S test returns a null probability of p = 0.998. 
Figure [7] shows the comparison of [Om] and optical lumi- 
nosity to the log 10 (L([On])/L([Om])) ratio for these RLQ 
and RQQ subsamples. Assuming that both optical and [Om] 
luminosity are good tracers of AGN power, it is clear that 
the log 10 (L([Oll])/L([Om])) ratio for both populations de- 
creases strongly with i ncreasing optical a nd [Om] luminos- 
ity, in agreement with iKim et alj (|2006T ). However, in the 
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Figure 7. log 1 Q(L([Oll])/L([Olll])) ratio versus optical luminos- 
ity (top) and [Om] luminosity (bottom). The red points represent 
the RLQs and the grey points the RQQs matched in [On] and op- 
tical luminosity. The solid lines with the same colours are the OLS 
bisector fits for the two samples. 
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case of RLQs the log 10 (L([Oii])/X([Oiii])) ratio is higher 
indicating the same result as found using only the optical 
luminosity matched samples. 

In Figure [8] we show the comparison of the [On] emis- 
sion line properties of the RLQs and the matched sam- 
ple of RQQs. Using the [Om] luminosity as an estimator 
of AGN luminosity, the [On] emission excess remains for 
the RLQs. We have applied a one-dimensional Kolmogorov- 
Smirnov test for each of the histograms in Figure [5] between 
the RLQs and RQQs. The K-S test returns a probability of 
p = 9.90 x 10~ 12 , p = 1.37 x 10" 6 and p = 5.32 x 10" 10 
and a test statistic of D = 0.28, D = 0.15 and D = 0.26 
for EW [0 ii], log 10 [L([Oll])] and log 10 [L([On])/L op t] respec- 
tively, thus the null hypothesis, that the three distributions 
are drawn from the same underlying distributions, can be 
ruled out. Therefore, the matched samples in [Om] lumi- 
nosity and redshift demonstrate that our results are ro- 
bust against correlations between radio luminosity and [Om] 
emission and the corresponding accretion rate derived from 
these quantities. 



4 DISCUSSION 

In this section we discuss the mechanisms that may con- 
tribute to the observed difference in [On] emission between 
the RLQ and RQQ populations. 

4.1 Mergers and triggering of activity 

The early claims that quasars have a high incidence of 
close companions led many studies to suggest a connec- 
tion between quasar activity, i nterac t ions and mergers (e.g. 
Urrutia et all I2OO8I ; IShankarl |2010| ; H amos Almeida et aJj 
20121 ). It has been suggested that major mergers of gas- 
rich galaxies are an efficient mechanism to feed supermassive 
black holes (SMBHs). This suggestion is supported by obser- 
vati ons of quasar hosts undergoing mergers and interactions 
(e.g. iKauffmann et all [20031 ; lUrrutia et I all 120081 ). Although 
we note that other deep field studies of the hosts of mod- 
erate luminosity AGN have found no significant evidence 
for high merger rates or interactions (|Sanchez et al.l [20041 ') 
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suggesting that secular processes may trigger these AGN 
(e.g. bars, slow cold gas accretion, minor interactions, su- 
pernova explosions: for reviews, see lKormendv fc Kennicuttl 
12004 lMartinill2003 ; |jogee!l2006l ). Mergers or inflows of cold 
gas by other mechanisms not only feed the black hole, but 
provide a new reservoir of available gas capable of producing 
a sudden burst of star formation. Based on th e scenario of a 
direc t link between mergers and quasars (e.g. I Sanders et al.l 
1988), strong interactions or major mergers between gas-rich 
galaxies fuel starburst activity and als o provide some of the 
gas for the black hole accretion event. [Bennert et al.l (|2008f ) 
suggest that most QSO host galaxies are involved in mergers 
with accompanying starbursts but usually several hundreds 
of Myr pass after the merger until the activity is triggered 
(see also lHopkins et al.ll2006al : IWild et al.ll20ldf ). This would 
explain the similarities in [On] emission for the RLQs and 
RQQs, at least at high luminosities. 

Major galaxy mergers are also expected to lead to black 
holes with a high value of spin while minor mergers tend to 
lead to low values of spi n. According to the spin paradigm 
l| Wilson fc Colbertlll995h the jet production and power are 
correlated with the angular momentum of a spinning black 
hole. The spin paradigm has often been used to explain the 
radio loudness dichotomy, implying that RLQs have black 
holes that spin more than RQQs. 

Based on the above theoretical background, mergers 
seem to be associated with the AGN activity and the black 
hole spin which is used to explain the radio loudness di- 
chotomy. However, can mergers on their own explain the 
[On] emission difference between RLQs and RQQs at low 
optical luminosity and the similarity in [On] emission at 
hig h optical luminosit y as found in Section 0? 

iDotti et al.l j2010l ) have proposed a model based on gas- 
rich mergers between nearly equal mass disc galaxies, 'wet 
mergers', that result mostly in radio-quiet AGN residing in 
elliptical or spiral galaxies with 'extra light'. The available 
gas and the interactions trigger a starburst to create the 'ex- 
tra light' at the center of the galaxy. On the other hand, ma- 
jor gas-poor mergers, 'dry mergers', are expected to lead to a 
radio- loud AGN hosted by an e lliptical galaxy. In agreement 
with this, iDunlop et all [(20031 ) have shown that radio-loud 
AGN are predominantly hosted by elliptical galaxies. How- 
ever, our results do not support this, because if RLQs are 
the result of major gas-poor merger then the available gas 
to trigger a starburst is not enough and the estimated [On] 
emission excess must be the result of a different mechanism. 
One possibility is that the high spin black hole produced 
by the merger powers the radio jets which compress the in- 
terstellar gas, triggering star-formation activity during the 
QSO phase. 



4.2 Alignment effect and positive feedback 

One important aspect of AGN is the role that they play in 
possible triggering of extens ive star formation i n a multi- 
phase intergalactic medium (|Silk fc Nu sser 2010). Star for- 
mation triggered b y expanding radio gala xy lobes, especially 
those of the FR II jFanaroff fc Rilevl 19741 ) type, may explain 
the alignment between large scale optical emission and the 
axes of their relativistic plasma je t s (e.g. iMcCarthvl 1 19871 : 



galaxies and RLQs are dr awn from the same parent popula- 
tion (e.g. lAnton ucci 1993j). It is believed that the difference 
in these two populations in terms of their observed proper- 
ties is the line of sight, which in the case of RLQs is closer 
to their emission axis jBarthellll989T ). Consequently, RLQs 
should also exhibit the alignment effect, albeit spatially un- 
resolved due to the proximity to our line-of-sight. 

There are three main physical interpretations for the 
alignment effect. The optical emission is proposed to be 
either: 1) dominated by light from young massive stars 
that formed when an expanding radio source (e.g. jets) 
collapse the dense gas providing add ed pressure (e.g. 



iBicknell et al] l2000l : iFragile et all |2004| ). 2) scattered ra 
diation from the AGN at the centre of the radio g alaxy 
l|di Serego Alighieri et all Il989l : iTadhunter et"ai1 Il992l ). or 
3) due to Inverse Co mpton scat tering of Cosmic Microwave 
Background (CMB) (|Dalvlll992l ). 

The first of these, where the radio source propagates 
into the interstel lar and intergala ctic medium can also pro- 
duce shocks (e.g. IBest et al.| [2000). Low-ionization emission 
lines could be produced by the ionizing photons in the cool- 
ing - shocked gas or, alter natively, by the warm clouds as 
they c ool b ehind the sh o ck. Solorzano-Inarrea fc Tadhunter] 
l|2003h and iTilak et all (|2005l ) have confirmed that the ra- 
dio lobes and hot spots are preferentially associated with 
lower ionization gas produced by precursor gas ionized by 
the shock being driven into t he cloud by the deflect ed radio 
jet. However, as discussed by 



ing lElmegreen fc Elmegreer 



Bicknell et all (|2000l ). follow- 



Dy tJickn 
iil (|l978l ) 



the shock ionized gas 



IBest et alj 1 19961 : lEales et all Il997l : ICroft et al.ll2006l ). Un- 



der simple orientation-based unified schemes, powerful radio 



also tends to go hand-in-hand with the occurrence of star- 
formation as gravitational instability occurs on timescales 
comfortably within the dynamical timescale of the jet-cloud 
interaction. Therefore, although we cannot rule out purely 
low-ionization shocked gas, it is unlikely that this is solely 
resp onsible for the excess [O n] emission. 

iBattve fc Brow ne (2009) compared the alignment effect 
for different radio-loud AGN at 1.4 GHz. They found that 
less-luminous RL AGN tend to align their radio emission 
axis with the minor axis of the starlight of the host. In con- 
trast, no significant preference of radio-optical alignment is 
found for more powerful RL AGN. This study appears to 
be at least consistent with finding that there is higher [On] 
equivalent width in RLQs with lower radio luminosity. 



4.3 Negative AGN feedback 

Semi-analytic models invoke AGN fe edback to truncate 
star formation in massive galaxies (e.g. ICroton et al. I l2006l : 
I Bower et al.l 120061 ). often termed negative feedback. This 
feedback could take several forms such as the photoioniz- 
ing and heating of the gas reservoirs in the host galaxy (e.g. 
iPawlik et al. 1 12003 ). or the exp ulsion of significant am ounts 
of gas through outbursts (e.g. iNesvadba et al.l [20061. The 
correlation between the formation of black holes and the star 
formation can be interpreted in two ways. In the first one, 
black hole accretion and star formation occur simultaneously 
because they are both fed from the same gas, brought to the 
center by gas-rich mergers and disc instabilities. Black hole 
accretion is then terminated when star formation has used 
up all the gas. In the second interpretation, star formation 
terminates when the black hole blows all the gas outside it s 
host galaxy (e.g. ISpringel et alll2005l : iHopkins et alll2006ah . 
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An interesting observational result is that most massive 
galaxies since redshift z ~ 2 appear to shut d own star for- 
mation earl ier than the less massive ones (e.g. iBundv et al.l 
l200rj. 2008 ). Thi s phenomenon, known as 'downsizing' 
( Cowie et al .1 19961 ). coincides with the drop in the space den- 
sity of quasars at recent epochs pointing to the possibility 
that black hole accretion and star formation are terminated 
together. 

4.4 A possible explanation 

The results of our analysis, and especially Figure|6] could tell 
a story about the star formation and [On] emission in RLQs 
and RQQs. The break at optical luminosity /black hole mass 
versus [On] equivalent width for the RLQs provides strong 
evidence for different evolution of star formation between 
the two quasar samples. It seems that when the black holes 
are less massive, presumably due to observing them ear- 
lier in their lifetimes, the two populations have significant 
differences in their [On] emission, but they have the same 
distribution when they have accumulated the bulk of their 
final mass. 

We propose that both of the populations start their ac- 
tivity due to a merger event or through an inflow of cold gas 
due to secular processes. As long as strong gravitational in- 
teractions play an important role in the production of super 
massive black holes and the quasar activity, this should be 
common for both radio loud and radio quasars. Depending 
on the pair of galaxies that were involved in the interaction 
and whether they can spin up a black hole, RLQs or RQQs 
could both be produced. For example, a high-mass/low-mass 
pair will produce a high mass black hole with low spin yield- 
ing a RQQ. On the other hand the merger of two high-mass 
galaxies would produce a high-mass, high-spin black hole 
and probably a RLQ. 

Physically, right after the halo merger, a torque is ex- 
erted in the host and gas is quickly funneled into the center 
of the newly-formed gravitational potential. In this phase, 
the black hole accretes its fuel very efficiently and the lumi- 
nosity of the quasar increases. In order to explain the [On] 
emission difference between the two populations, we need 
a mechanism which can trigger star formation in the case 
of RLQs. The obvious mechanism in this case are the jets. 
As the jets pass through the emission line region, the shock 
will provide additional ionizing photons, increasing the [On] 
emission line. Radio jets are able to induce the formation of 
massive knots of bright young stars, which will disperse and 
fade over the lifetime of the quasar. This mechanism results 
in the tight alignment of the bright optical-UV continuum 
emission along the radio jet. 

As the black hole mass increases, the jets are no longer 
able to keep the gas away from the black hole. A galaxy with 
a large bulge has been created and the new stars can feed 
the black hole. At this point the negative AGN feedback is 
stronger than the positive, the accretion rate and the black 
hole mass increase and as a consequence the luminosity of 
the quasar. The [On] emission in RLQs decreases relative 
to the RQQs, and the RLQs and RQQs would have similar 
[On] emission. For both of the populations the accretion can 
not last indefinitely, and it is quenched when the quasar is 
powerful enough to expel this gas or at least to balance its 
infall, thus inhibiting further growth. 
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Figure 9. Comparison of the EW([On]) distribution for the sam- 
ple of RLQs (red) and the matched sample of RQQs (black). The 
dashed lines shows the mean values for each population. The top 
panel shows the distribution for sources with log 10 (L op t/W) > 
38.6 while the bottom for sources with \og 10 (L opt /W) ^ 38.6. 



In order to support the above explanation, in Fig- 
ure we compare the [On] emission for the two popula- 
tions for optical luminosities log 10 (I/ op t/W) > 38.6 (top 
panel) and log 10 (L opt /W) < 38.6 (bottom panel). The 
difference is obvious for the [On] equivalent width. For 
log 10 (Lopt /W) > 38.6 the mean [On] equivalent widths for 
the two populations are < EW([Oii])r L >= 4.65 ± 0.05 A 
and < EW([Oii])rq >= 4.83 ± 0.04 A for RLQs and 
RQQs respectively. In contrast, for log 10 (L opt /W) ^ 38.6 
the mean values are < EW([Oii])rl >= 8.74 ± 0.08 A and 
< EW([Oii])rq >= 4.67 ± 0.05 A. The one-dimensional 
K-S test returns for log 10 (L op t/W) > 38.6 a probabil- 
ity of p = 0.24 and a test statistic of D = 0.08 while 
for log 10 (L opt /W) sj 38.6 returns p = 1.04 x 10~ 16 and 
D — 0.18. For high optical luminosities the two populations 
are drawn from indistinguishable distributions and have the 
almost the same mean [On] equivalent width, whereas for 
low optical luminosities RLQs have approximately twice the 
[On] emission of RQQs. As we can see, in the case of RQQs 
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there is no significant difference in [On] emission for the two 
optical luminosity bins. 



5 CONCLUSIONS 

In this paper we have studied the relationship between the 
[On] emission line properties of matched samples of RLQs 
and RQQs. The main result of our study is that RLQs have 
higher [On] emission than RQQs. We find that the [On] 
equivalent width is higher for RLQs at low optical luminosi- 
ties. As the optical luminosity and also the black hole mass 
increase, the [On] equivalent width for RLQs decreases while 
for RQQs it remains almost the same. Above an optical lu- 
minosity limit of log 10 (L op t /W) = 38.6, the [On] equivalent 
widths for the two samples are indistinguishable. Below this 
limit the L opt — EW([On]) correlation changes for RLQs, 
while for RQQs it remains the same. 

Based on these results, we believe that the main reason 
for this difference is the presence of powerful radio jets in 
RLQs that pass through the gas and delay the in-fall of the 
gas into the central supermassive black hole. At the same 
time they ionize the gas in the ISM (interstellar medium) 
and are able to induce star formation. When the black hole 
mass increases, along with the luminosity, due to the accre- 
tion of matter, the jet power is no longer enough to keep 
the gas away from the black hole and the star formation 
decreases. At these high optical luminosities the two popu- 
lations tend to have the same [On] emission. However, we 
note that to obtain a full picture of the process occurring 
then a measurement of the environmental dens i ty would also 
be advanta geous (e.g. iKauffmann et aL I l2008l ; iFalder et al.l 
l20ld , l201lD . 

Even if [On] emission is generally used as a star forma- 
tion tracer, it is not one of the strongest indicators and our 
results may be biased due to the AG N emission especially at 
high optical luminosities (|Holl2005h . although we emphasize 
that our two populations have the same optical luminosity 
distribution. However, the fact that we show that the [On] 
emission excess for the radio-loud population is present using 
either the [Oin] luminosity or the optical luminosity as an 
AGN luminosity tracer, reinforces our suggestion that this 
excess is the result of higher star-formation activity in RLQs. 
A caveat to this is that some, possibly significant fraction 
of the [On] emission could emanate from shocked gas within 
the radio lobes, although we note that such emission usually 
goes hand-in-hand with star formation. Our results therefore 
open the way for investigations in this direction using more 
secure star formation tracers such as far-infrared emission. 
This will be possible with surveys such as Herschel- ATLAS 
which cover a wide enough area to include a large number 
of the rarer RLQs. 
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